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We phenomenologically developed a propagation model of high energy galactic cosmic 
rays. We derived the analytical solutions by adopting the semi-empirical diffusion equa- 
tion, proposed by Berezinskii et ai. (1990) and the diffusion tensor proposed by Ptuskin 
et aZ.(1993). This model takes into account both the symmetric diffusion and the anti- 
symmetric diffusion due to the particle Hall drift. Our solutions are an extension of the 
model developed by Ptuskin et al. (1993) to a two-dimensional two-layer (galactic disk 
and halo) model, and they coincide completely with the solution derived by Berezinskii 
et al. (1990) in the absence of antisymmetric diffusion due to Hall drift. We showed that 
this relatively simple toy model can be used to explain the variation in the exponent of 
the cosmic ray energy spectrum, 7, around the knee E 10-*^^ eV. 
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1. Introduction 

Cosmic ray propagation is one of the most important and interesting subjects in 
astrophysics and high energy particle physics. It is believed that the observed cos- 
mic ray data such as the cosmic ray energy spectrum includes information about 
the space through which the cosmic rays pass. In fact, it is possible to evaluate the 
thickness of matter, and it is also thought that details about the galactic magnetic 
field can be extracted from the observed data because the cosmic rays experience 
frequent collision and scattering with both the interstellar gas and the galactic 
magnetic field during their propagation in the galaxy. Thus, a reliable cosmic ray 
propagation model may enable us to obtain further knowledge about the galactic 
structure. Thus far, several cosmic ray propagation models have been proposed and 



1 



December 8, 2010 1:39 WSPC/INSTRUCTION FILE arakida-ijmpa 



2 H. Arakida and S. Kuramata 

discussed (see and references therein). Further, some researchers have success- 
fuUy derived analytical solutions using the diffusion equation SSEI^. 

The spectrum of cosmic rays shows one of the most distinctive features, known as 
"knee" , around energy E w 10^^ eV at which the exponent of the energy spectrum, 
7, changes from 2.6 - 2.7 for IQio <E< 10^^ eV to 3 - 3.1 for lO^^ <E< 10^® eV. 
Currently, it is not clear why the exponent changes around E « lO^'^ eV. Thus far, 
several models have been proposed to explain this spectral property: a shock wave 
acceleration model based on the acceleration of cosmic ray particles by the shock 
wave front ^ ^"^j a diffusive propagati on model based on the leakage and 

diffusive propagation of cosmic rays in the galaxy ^ '^ l -^^^ l ^ l ^^ l ^^ l ^^ l -*^^ ^^. an interac- 
tion model based on th e interaction of cosmic rays with the background particles in 
the galaxy | 21 | 22|23 |24j^ a rea ction model based on the reaction of cosmic rays with 
the atmosphere of Earth l^^^^fi,^ g|-^^ Among these models, the shock wave accelera- 
tion model seems to be widely accepted as the explanation for the knee. However, 
the diffusive propagation model can provide the exposition for th e first kriee and 
the second knee and for the observed compositions and anisotropics ^^2^11^1^. The 
diffusive propagation model is characterized by introducing the particle Hall drift 
effect; thus, it is a theoretically simple model. 

In this study, we derived solutions for the diffusive propagation model of cosmic 
rays and confirmed the validity of this model. We adopted the semi-empirical dif- 
fusion equation introduced in ^ and the diffusion tensor described in El. Then, we 
extended the propagation model in 1^ to a two-dimensional two-layer (comprising 
the galactic disk and the halo) model in cylindrical coordinates. 

This paper is organized as follows: In Section [2l we derive the analytical solutions 
for the diffusion equation of cosmic rays. In Section [Sj we qualitatively show that 
our model can explain the spectral feature of the observed cosmic rays, namely, the 
exponential variation around the knee. In Section |4l we present the conclusions of 
our study. 



2. Solutions for Diffusion Equation 
2.1. Galactic Structure 

Fig-fflshows a schematic diagram of our galactic model. We assume that our galaxy 
has a cylindrical structure and it is divided into two parts; the galactic disk (D) and 
the halo (H). The signs + and — indicate the direction along the z-axis. In Fig.[Tl 
R is the radius of the galaxy and 2hg and 2h are the height of the disk region and 
the galaxy, respectively. We suppose that the source of cosmic rays is distributed 
uniformly within the galactic disk area only. The source region of the i-th particle 
is given by radius Ri {0 < Ri < R) in the shaded area in Fig. [TJ 
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Fig. 1. Schematic Diagram of Galactic Structure. 



2.2. Basic Equation, Diffusion Tensor, and Hall Drift 

We consider the following transport equation 

- Va{Dafi{z)V0)N,{E, r, z) + n{z)va,N,{E, r, z) = Q,{E, r, z), (1) 

where Ni{E,r, z) is the number density of the j-th cosmic ray component per unit 
time, unit volume, and unit energy; Z3q^ is the diffusion tensor (indexes a and 
/3 denote the coordinates); v is the velocity of the cosmic ray; ai is the inelastic 
scattering cross-section between the i-th cosmic ray particle and the interstellar gas; 
n{z) is the density of the interstellar gas; and Qi{E, r, z) represents the distribution 
of the source. 

We adopt the diffusion tensor given in ^ 

Dap = - D±_)babfj + D±5al3 + DAeapjb^, (2) 

where ba is the unit vector of the magnetic field, -D|| denotes the diffusion along 
the magnetic field line, D± denotes the diffusion perpendicular to the magnetic 
field line, Da is diffusion due to the antisymmetric drift of the particle, 6a/3 is 
the Kronecker's delta symbol, and ea/s-y is the complete antisymmetric Levi-Civita 
tensor. 

The observed data show that the magnetic field in the galaxy has a spiral struc- 

F97I 

ture, and it is regarded as an approximate toroidal structure'^. Thus, we assume 
the unit vector of the magnetic field as 

6, = 63 = 0, b,p = 1. (3) 

It should be noted that the magnetic field line in the galactic disk changes the 
sign almost every 3 kpc along the radial direction Nevertheless, for the sake of 
simplicity, we assume that the global magnetic field points in the same direction 
anywhere in the galactic disk. On the other hand, the structure of the magnetic 
field in the galactic halo is not well known; so we suppose that the magnetic field 
in the halo also has a toroidal structure given by 

= = &f ) = = 0, bl^^ = b'-P = 1, (4) 
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where (D) and (H) denote the disk and the halo, respectively. From Eq. (jH), the 
diffusion terms in ([ij becomes 

Id d d d d did d 

- V.P^.V^) = —UrD^)^ ^D^^ ^Da-^ + -UtDa)^. (5) 
r dr dr dz dz dz dr r dr dz 



These coefficients are determined in0 



as 



D\\ = y , D^= gA — , Da = —■ (6) 

Here I is the mean free path that is characterized by the turbulent exponent m, A 
denotes the relative value of a random magnetic field within the characteristic scale 
L, and g is a not well-determined parameter (see^. From the gyro-radius 

c B xp 

^^^0 = ^-^' (7) 



we obtain the drift velocity of the particles 



pcv Bj^ 

Vd-i = ~Vy{DAeafiyba) = T^^ap-r^—-^- (8) 



Then, we substitute the coefficients of the diffusion tensor given in 



T 

D±{r,z) = constant, £'^(r, z) — Dao{z) — , Dao{z) = constant. (9) 

K 

Then the diffusion term reduces to 

V7 m T7 ^ D±d d d^ 2Dao d 

r dr dr dz-' R dz 

From Eq. ([9]), the drift velocity Vd is written as 

Vor = ~-^Da = 0, Vd^ = --^{rDA) = = constant. (11) 

dz r dr R 

Thus, the drift effect appears only along the z-direction in this case. 



2.3. Boundary and Continuity Conditions 

Let us summarize the boundary and continuity conditions. First, at the edge of the 
galaxy, we assume that the cosmic rays leak from the galaxy 

Ni{r = R,z) = N,{r, z = ±h) = 0. (12) 



^ It is suggested that the turbulent exponent in the magnetic fiel d in th e galactic disli is practically 
characterized by the power law with exponent m = 0.3 ± 0.3 I^^ I ^ H For instance, m = 1/3 is 
related to the Kolmogolov turbulence spectrum; m = 1/2, t urb ulence spectrum of Kraichinan 
hydromagnetic one ED and m = 0, Bykov-Toptygin spectrum 1221. In the case of the halo, we do 
not have sufficient knowledge about the magnetic turbulence; therefore, we assume that the same 
power law holds 
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This condition seems to be appropriate in terms of the cosmic ray L/M ratio. 
Second, at the boundaries of the disk and the halo, z = ±hg, both Nj;'^^ and ivl^'' 
must be continuously connected. Then, 

and 

dz dz 

2.4. Solution for Transport Equation 

Because Eq. ([T]) has the same form as Eqs. (3.10) or (3.11) in 131^ we follow the same 
approach. To obtain the solution for Eq. ([Ij, wc first obtain the Green's function 
$(r, z; ro), which satisfies 

- V^{D^p{z)Vp)^{r, z- ro) + n(z)w,$(r, z; tq) = Qihg - \z\y^^—^, (15) 

where ro represents the position of the cosmic ray source and d(h,g — |z|) is the step 
function. $(r, z;ro) is related to \y^. 

N, ^ ATThgfj.E-''^'' xiroMr,z;ro)rodro, (16) 
Jo 

where gi is the constant characterizing the j-th particle, x(''o) represents the radial 
distribution of the source of cosmic rays, and we presume the energy spectrum at the 
source obeys the power law E~^^° for the i-th particle. We also suppose the radial 
component of $(r, z;ro) is expanded by using the zeroth-order Bessel function 

oo 

$(r,z;ro) =^a,(z;ro)Jo (^) . (17) 

s — 1 

In this expression, Vs is the s-th root of the equation Jo(i^s) — 0. Hence, Eq. (jlSp 
becomes 

2^ 1 —d-/d-r - ^^("te + -^d;\ + "("^"^7 «4^;ro) Jo j 

s— 1 ^ 

= eihg-\z\)'±-^. (18) 

47r/igro 

From the orthogonality of the Bessel function, the delta function is expressed as 

*-™)^£p^i*(^)-'»(¥)' 

where Ji is the first-order Bessel function. Noting that 
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Eq. is written as 



2 92 2Dao{z) d 



.R. 
9{h,~\z\) 



R dz 



{z)vcji \ as{z; ro)Jo (-^) 



Thus, the problem reduces to solving the ordinary differential equation 



R 
(21) 



2Dao{z) d 
dz^ ^ RD±{z) 



where 



C 



as{z;ro) - 
e{h,-\z\) 



!^s\2 n{z)va. 



27r/igi?2Z?(2)^|Ji(i/,)|2 
Daq, and n(z) are chosen as 



D±{z) 



a,(z;ro) = -C, (22) 



(23) 



D^{z) = 



Da 



Dao{z) = 



Dao 

dAO 



n{z) 



Ug \z\ < /ig(disk), 
\z\ > hg{halo). 



In the galactic halo, C = because of the step function. Therefore, 



d'ai"^ 



2dAo da 



(H) 









d± \ 



(H) 



dz"^ Rd±_ dz 
and the general solution has the form, 

ai^+)(z;ro) = Ae"'^ +Be'5'^ z > 0, 



(24) 



(25) 



(26) 
(27) 



where A, A' , B and B' are the integration constants and at and /3i are the solutions 
of the following equation that is associated with Eq. 





Hhvai 




di_ \ 



Rdi 



Then the following two solutions are obtained 

dAO 



= 0. 



^3^ = 



Rd± 

dAO 

Rdl 



f dAa\ , (^s\^ , nhVUi 
-^[-Rdl) +[r) 

, // dAo\^ , f'^s\^ , UhVOi 

^n-Rdi) +[r) 



(28) 

(29) 
(30) 



From the boundary condition in Eq. (|12p . Eqs. (|26l) and (l?7l) are rewritten as 

ai"+\z;ro) = B[-~exp{{(3,-a^)h}e°'^' +6^^'] z>0 (31) 



\z;ro) - S'[-exp{-(/3» -a,)/i}e"'^ +6'^''=] z<0. 



(32) 
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In the galactic disk, the general solution takes the form 

ai°+^ (z; ro) = Pe"''' + Qe^'' + C z>0, 
a(^-) (2; ro) = P'e^'" + Q'e**" + C z < 0, 

where P,P',Q and Q' are also integration constants and 

C' = 



C 



"fi and Si are the solutions to an equation similar to Eq. ([28l 
Dao 



Dao 



Dao Y /i^y 
RDj_ J \r) 



UgVCTi 

D± ' 



We suppose the cosmic ray source distributes uniformly; then 
x(ro) = 9{R, - ro), < R, < R. 

Furthermore, 



R^ 



0{Ri - ro)ro Jo (-^) dro 



RRi 



R 



(33) 
(34) 



(35) 

(36) 
(37) 

(38) 

(39) 



and after straightforward but bit tedious calculations, we obtain the number density 
of the i-th cosmic ray particles 



E 



t|Ji(^.)|V.Z?^i?[(t)^ + ^ 
1 



g{Si-li)hg _ g-(<5i-7i)'»g 



{lii-a,){h-h,) 



[(7i - a,) exp{(/?, - a,)h}e'^^^^ - (7^ - P^)eP^^^] 

X ((*» - + ml^ + 5i) - 2yA ~ (a? + fi^^)] 

(a, - - A)e^''~^'^'*'' - (A - - a.)e-'''-^""''' 

(A - - aOe^''-^'^"' - {a, - 5i){lr - A)e"^''~^'^'« 

X { [(a, + /3,)7. " + A')] (e-^-'^' - e^^''^^ 

[a,(7, - /3,)e''''*« - ^(7, - a.)e-T'''»] e-(^'-"')(''-''«)}e^''*« - 7i)eT'''. 
x[-exp{(A-a,)Me"'' + e^'1, 



(40) 
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g 2giE-'>->R,.h (^) Jo (to 



1 



[(7, - ai) exp{(/3, - a^)h}e'^^''^ - (7^ - P^)ef^^<'] 
X ([(ai - 5i) exp{(A - (Xi)h)e^^''o - (A - k)e^''' 
X { [{ai + + Si) - 2jA - {a^i + /Sf)] 

- [{ai - 6i)i^i - A)e(*^-^*)'^^ - {Pi - Si){^i - a,)e-(**-^*)''='' 

- [(ft - Si}{y, - a,)eS'^-'^^^^^ - {ai - (5,)(7. - ft)e-(**-^*^'^''' 

X { {{ai + A)7i - (a? + m\ (e^-'^" ' 

- [ft(7i - ai)e^'^'> - ai{-n - ft)e-^-''«] e(ft-".)(''-'>«) 

+ai exp{(A - ai)/i}e«*''=' - Ae^*''=')eT*^ 



-(ft-ai)(fc-/tg) 



+ 1 [(a, + ft)(7« + <5.) - 27.<5. - (aj + )] 

- \ai - ^i)(7i - A)e(''-^''"'' - {Pi - 5i){j, - a,)e-^''-'^^''^ 

- [{Pi - 5i){-/i - aOe^^^-^')'^" - {ai - 5i){^i - pi)e-^''-^'^^^ 
X { [(ai + Pihi - (a? + Pf)] (e^^''^ - e"^*''") 

- [Pi{7i - ai)e^'^^ - ai(7i - Pi)e-^'^^] e^0i-'-')ih-hg) 

- [ai{li - ft)e^-'^« - A(7i - ai)e-T-''«] e-^P^-<-^)(h-hg)ys.. ^ 



-(13, 



-ai){h-hg)^ ^ 



(41) 
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E 



(H-) 



{E,r,z) 
25,i?-^»i?,Ji(i^)Jo(tr) 



1 



[(7i 



{0i-ai){h-hg) 



{I3i-<xi){h- 



- [(A - 5i){'yi - ai)e^^'-^'^^^ - (a^ - J,) (7. - P^)e-^^'-^'^''■ 



('5.-7.)fts 



x[- exp{-(ft - a,)/i}e"- 



(42) 
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nI"~\e, r, z) 



1 



[(7, - a,) exp{-(ft - aO/ije-"-'*-'' - (7^ - P^)e'P^^^] 
X { [(a, + A)(7^ + Si) - - (A + [e(«'-^')''« - e 

X { [(a. + A)7. - (a- + /3- )] (e'^''^^ - e-^-''^) 

+ { [(a, + ft)(7. + Si) - 27.<5. - {Oil + /?f )] [e^^'"^')''- - 

(a, - (5.)(7. - A)e(^'-^')'*« - (A - 5,)(7^ " aOe"'''-^-"''' 

x{[(a, + A)7. - (a- (e-^'"' " e^^-'^O 



^{P.-a.){h~h,) 



-(/3.-a.)( 



(.5. -7.) 'is 



-(ft 



-aOCi-'i.,)! ^ 



(43) 



We note that these solution reduce to those in ^ in the absence of antisymmetric 
diffusion due to particle drift. 

3. Exponential Variation of Cosmic Ray Spectrum 

In this section, we discuss the features of the derived solutions, especially focusing 
on the question whether the diffusion model can reproduce the exponential variation 
of the cosmic ray spectrum. Here, we present qualitative examples and then omit 
the scales in the following figures. 

In accordance with we assume the diffusion coefficients have the following 
energy dependence 

D^ = Dd±E'^^, d^^dH^E"^", Dao = DdaoE, dAo = dHAQE, (44) 
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Fig. 2. Energy (E) Dependence of Intensity li. 

where the turbulence exponents in the disk and halo are denoted by m/j and mn, 
respectively, and for the energy spectrum at the source, we choose 7,0 = 2.2. 

Fig. [2] shows the energy {E) dependence of the intensity of the model by ^ 
(cited as Ptuskin) and out three solutions: (1) Du^ = dHi_ and Ddaq — dnAo 
(shown as = dn)', (2) the diffusion coefficient of the halo is larger than that 
of the disk, D]j±^ < dHi_ and -D_dj^ < dni, (shown as < dn)', and (3) the 
diffusion coefficient of the disk is larger than that of the halo, -D_dj_ > dHi_ and 
Dd± > dH± (shown as Do > dn)- In the case of our solution, the exponents of 
magnetic turbulence in the disk and the halo are equivalent, mu = mn- In Fig. 
[21 we multiplied E'^'^ with li to emphasize the exponent variation. The three cases 
— Du = dn, Du < dn and Du > dn — can reproduce the variation in the 
cosmic ray spectrum in the same way, as the model developed by Ptuskin. Fig. |3] 
shows the m dependence of our solutions: (1) the exponents of the disk and the halo 
are equivalent, ttid = mn', (2) the exponent of the halo is larger than that of the 
disk, mo < ruH', and (3) the exponent of the disk is larger than that of the halo, 
iTLD > tuh- In all cases, we fixed Db±^ — dH± and Duao = duAo- We found that all 
the solutions exhibit a similar trend, as shown in [2j Further, we may say that the 
difference between the two exponents around the knee, A7 = 7_e<_e^„^^ — 7_Ek„„„<_E, 
is larger value for mo < rnn than that for mr) > mn- 

4. Conclusions 

We phenomenologically proposed a propagation model of galactic cosmic rays based 
on the semi-empirical diffusion equation developed by Berezinskii et al. (1990) and 
the diffusion tensor introduced by Ptuskin et al. (1993). This model takes into 
account both the symmetric diffusion and the antisymmetric diffusion due to the 
particle Hall drift. The derived solutions are an extension of the model developed by 
Ptuskin et al. (1993) to a two-dimensional two-layer (galactic disk and halo) model. 
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Fig. 3. E Dependence of Intensity li . 

and they coincide completely with the solutions derived by Berezinskii et al. (1990) 
in the absence of antisymmetric diffusion due to particle drift. We shown that this 
relatively simple model can be used to explain the variation in the exponent of the 
cosmic ray energy spectrum, 7, around the knee E w IQ-'^^ eV. 

In this paper, we showed that although the diffusive cosmic ray propagation 
model can be used to explain the observed cosmic ray spectrum, especially the ex- 
ponential variation around the knee, our mode;! is actually a more simple toy model 
based on assumptions such as the cylindrical structure of the galaxy and simplifica- 
tion of magnetic field. To further test the validity of the diffusion model, we must 
conduct numerical simulations under more realistic situations. This is a difficult 
task; nonetheless, it may help us to gain a deeper imdcrstanding of astroparticle 
physics and the galactic structure through which the cosmic rays pass. 
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